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Abstract 

We present a treatment of cold hydrogen-antihydrogen collisions based on the asymptotic prop- 
erties of atom-antiatom interactions. We derive general formulas for the elastic and inelastic cross 
sections and for the scattering lengths and analyze their sensitivity to the parameters characterizing 
Q_i, the inelasticity of the collision process. Given the inelasticity, we obtain bounds for the complex 



O 



scattering length. We investigate the influence of strong nuclear forces and the isotope effects in 
HH and HD collisions and demonstrate enhancement of these effects due to the presence of the 
near-threshold narrow HH (HD) states. The values of the elastic and inelastic cross-sections with 



simultaneous account of rearrangement and strong forces are presented. General expressions for 
the (complex) energies of the near-threshold HH states are obtained. 
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I. INTRODUCTION 

A very important property of ultra-cold atom-antiatom collisions is the separation of 
scales, characteristic for inelastic transitions on one hand, and van der Waals interaction on 
the other hand. Indeed, the analysis of the dynamics of Protonium (Pn) and Positronium 
(Ps) formation [l|, [2|, y, |4j, ly, la, L3, la, l9( , which is the dominant inelastic process in H — H cold 
collisions, shows that such transitions happen mainly at internuclear separation distances 
less than R s , which is on the order of the Hydrogen Bohr radius (re). As shown by previous 
investigations [6j this distance is optimal for the overlap of the initial channel (H — H) wave 
function and the final state (Pn — Ps) wave-functions. 

At separation distances larger than FL the H — H interaction is well approximated by 
the one-channel adiabatic potential V a d |Sl, llO] • Therefore the wave-function for the relative 



motion of the H — H could be obtained from the solution of the one-channel Schrodinger 
equation with the potential V a d(r), together with a proper boundary condition imposed at 
the internuclear distance r ~ R s . Importantly, such a boundary condition (which ideally 
should be derived from the solution of the four-body problem at internuclear distances 
r < R s ) turns out to be energy independent for sufficiently small energies of colliding H 
and H . Indeed, as follows from the kinematics of Pn and Ps formation 7{] the energy of 
relative motion of Pn and Ps is e = E + M/(2N 2 ) — m + m/(4n 2 ), where E is the H — H 
cm. collision energy, M is the reduced mass of Pn, m is the reduced mass of Hydrogen, 
TV is the principal quantum number of Pn, n is the principal quantum number of Ps. As 
one can see the minimum of e is e m = 0.046 a.u. (E = 0, N = 24 and n = 1). Therefore it 
is reasonable to expect that for H—H collision energies E <C e m the Pn-Ps wave-function 
would weakly depend on the collisional energy E. In the following we will be interested in 
cold collisions at energies E < 10 -5 a.u. where the S-wave contribution is dominant and the 
above condition is easily fulfilled. Thus the influence of rearrangement and annihilation on 
the wave-function in the HH channel can be described by the complex boundary condition 
for the logarithmic derivative of the wave- function at R s , which is energy independent in 
the energy domain of interest. This boundary condition can be conveniently expressed in 
terms of the short-range complex phase-shift 5 = 5r + iSi, which is connected to the above 



mentioned boundary condition via the relation 

^ = P (R s )cot(5) (1) 

where p(R s ) = \f2MV a d(R s ) is a classical local momentum given at the distance R s . 

In our study we will be interested in the near-threshold properties of H — H scattering. 
Such properties are critically dependent on the existence of the near-threshold singularities of 
the ^-matrix, namely the weakly bound states or resonances. It is known that the density of 
the near-threshold levels is determined by the large distance behavior of adiabatic potential 
Vad{r) UJ- In our case this potential at large distances turns into the homogeneous van der 



Waals potential —Cq/t 6 . The strength of this potential (Cg) determines the characteristic 



distance R v dw = \/2MC§ ~ 10.5 a.u. which, as will be shown later, gives the order of 
magnitude of the scattering length for the H H system. At the same time the variation 
of the boundary condition (related to the short range phase-shift S, see eq. [1]) will only 
weakly affect the density of the near-threshold levels; it results mainly in shifting of the 
whole spectrum of such states. It may happen that for certain value of 5 there is a state 
(resonance) very close to the threshold. This results in a strong enhancement of the cross- 
sections. 

In this paper we will study the evolution of the scattering observables as a function of 
the short-range complex phase-shift S, thought of as a free variable. Though the "true" 
value of 5 could be in principle obtained from the solution of the four-body problem, it 
is very instructive to get a general dependence of the scattering amplitude on 5. Indeed, 
the existing model calculations of the low energy H — H scattering are all restricted by 
approximations in which the effect of certain decay channels are neglected. In terms of the 
here developed approach any account of additional channels, or generally any improvement 
of the optical potential, result in a change of the short range phase-shift 5. We present the 
universal dependence of the scattering length on the complex short-range phase-shift 5 in a 
closed form and show that such dependence manifests resonant behavior, connected to the 
existence of narrow near-threshold HH states. We analyze this behavior in terms of the 
near-threshold poles of the S'-matrix. Such resonance phenomena can significantly enhance 
the influence of "small" physical effects (or their neglect in approximate treatments) and 
make inapplicable naive perturbation approach. In particular, we study the influence of 
strong forces in the non-relativistic HH Hamiltonian, including the shift and splitting of the 



quasi-bound near-threshold HH states, as well as the isotope effect in D — H scattering. 

II. LOW ENERGY S-STATE H-H SCATTERING 

We start with the derivation of the analytical expression for the scattering length for H — 
H collisions. This derivation is based on the matching of the WKB form of the wave-function 
and the analytical zero-energy solution of the Schrodinger equation with the homogeneous 
(—Ce/r 6 ) potential in the asymptotic re gion . Such an approach was successfully used for 



the study of cold atomic collisions 
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131 ] . Our treatment will be extended to the atom- 
antiatom case, where account of inelastic transitions is important. 

Further analysis is based on the fact, that the WKB approximation is applicable for the 
description of the wave-function in between the short-range distance and the asymptotic 
domain R s < r <C R v dw'- 

i r 

$(r) ~ ^^=sin((5+ / p(r')dr') (2) 

where p(r) = \/2MV a d(r) is the classical local momentum in the adiabatic potential V a d(r), 
while 5 = 5r + iSi is the complex short range phase-shift, whose value should be fixed by 
matching with the solution of the four-body problem at r = R s - 

It is important that the distance (Rh) above which the adiabatic potential V a d(r) is well 
approximated by the homogeneous —Cq/t 6 potential is within the distance of the validity of 
WKB approximation, i.e. Rh < R v dw- The rough estimation of Rh can be obtained from the 
condition that for the distances above Rh the term C$/r 6 should dominate in the multipole 
expansion of V a d(r). This gives Rh = ^Cq/Cq ~ 5 a.u.. Thus there exists a matching region 
Rh < t <^i R v dw-, where both the WKB approximation and the zero-energy solution of the 
Schrodinger equation with homogeneous (—Cq/t 6 ) potential are valid approximations of the 
HH wave-function. 

In the vicinity of Rh, where the adiabatic potential is well reproduced by a homogeneous 
potential, the classical momentum is p(r) = ^2MCq/t 6 and one can get the explicit r 
dependence of the wave- function introduced in eq. (J2J): 

»( r )„ _^L sin( , + SJ - ±S (3) 



where we have introduced the semiclassical phase Q defined as 



Q = p(r)dr. (4) 

JRs 

The value of fl calculated using the adiabatic potential V a d of the H — H system from the 
ref. [14| and R s = 1 a.u. turns out to be fi = 19.383. 

The wave function given in eq. ([3]) can be matched in the vicinity of Rh with the exactly 
known zero-energy wave function <3>o( r ) m a homogeneous potential 1/r 6 [15], which is a 
linear combination of the form: 



*o(r) ~ V~r (j 1/4 (^p5) - CY 1/4 (- 



2r 2 

The coefficient C has to be determined from the matching procedure at R ~ Rh which 
results in the following wave-function at distances r > Rh'. 



*(r > R h ) ~ V~r (-M^P> ~ tan(| + (1 + S)Yy,C^)) . (5) 

Using the Taylor expansion of the Bessel functions for small argument and taking into 
account that the scattering length appears in the asymptotic form of the wave function 
through 

$(r — ► oo) ~ 1 — r/a (6) 

we obtain for the scattering length 

a = a (l + cot(| + Q + 5)) (7) 

where do is 

a ° = ^V5iwij" 4 '" a ' u " (8) 

As one can see R vc lw gives the characteristic size scale for the H — H scattering length. The 
argument of cotangent in ([7|) has a simple meaning. Indeed, Q is the semiclassical phase 
accumulated in the region r > R s , 5 is the short range phase-shift accumulated at r < R s 
and the term (7r/8) is the quantum correction to the semiclassical phase from the asymptotic 
van der Waals tail (— C 6 /r 6 ). 

The accuracy of the above expression is limited by the accuracy of the WKB approxi- 
mation in the range R s < r < R v aw and on possible discrepancies between the adiabatic 



potential V a d and the van der Waals term —C^/r 6 at the upper limit of that range. The 
comparison of (J7|) with the numerical calculation of the Schrodinger equation with the adia- 
batic potential V a d{r) and boundary condition (pQ) shows that formula ([7j) has the accuracy 
of 10 % in the wide range of 5. 

One can extend the above analysis to the calculation of the effective range. According to 
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17l | the effective range is: 



ao [ r(l/4) 
Te 3 [r(3/4)_ 

Upon substitution of a from eq. (JTj) one gets 



,o 



2^ + 2^ 



«o 



r(i/4) 



(cot(7r/8 + Cl + 5)) 2 + 1 



(9) 



e 3 [r(3/4)J (cot(7r/8 + fi + (5) + l) 2 ' 

For sufficiently low energy, i.e. when the momentum of the incident atom satisfies k\a\ <C 

1 the elastic (a e i) and inelastic (ai n ) cross-sections are determined by the scattering length: 
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;i + Recot(- + tt + S)) 2 + (Imcot(- + Q + S)f 



A7f\lma\/k = 47r^Imcot(- + Q + 5). 
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These values can be corrected to account for the second order terms in momentum k: 

Aii\a\ 2 
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1 - 2Hma+ k 2 (\a\ 2 - Re(r e a)) 
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k l-2k!ma + k 2 (\a\ 2 -Re{r e a))' 
Since the formation of Protonium and Positronium in the rearrangement collisions ends in 
complete annihilation, the total inelastic cross-section can be identified with the cross-section 
for annihilation. 

We will now turn to the study of the general dependence of the scattering lengths and 
cross-sections (elastic and inelastic) on the variation of the real (5r) and imaginary (5j) 
parts of the short range phase-shift S. 

Let us treat two important limiting cases. The first, which we refer to as the "weak 
absorption limit", is characterized by small imaginary phase-shift 5j <C 1 which means that 
the probability of transitions to the inelastic channels is small. In that case we get for the 
scattering length a from eq. (JTj) 



a f 1 + cot(- + CI + 8 R )J - ia Q 5i/ sm 2 (- + CI + 5 R ) 



(14) 



2 sin 2 (f + n + s R ) 

a eZ = 4vr|a| = 27ra . » 15 

S1I1 Z (| + il + R ) 

c 

cr in = 47r| Imal/fc = 47T- ^- -. (16) 

1 W ksm 2 (l + n + 6 R ) K ' 

The scattering cross sections (eqs. [ToTfToT) show oscillating character as a function of 5 R . 
As will become apparent later such oscillations originate from the existence of a spectrum 
of long-lived near-threshold states of the H H system. 

For the application to collisional cooling we are interested in the competition between the 
rate of elastic scattering and the rate of annihilation. The ratio cr e i/ai n reaches its maximum 

cr e A a k 



I ->, ( 17 ) 

(Jin / /Or 

m / max 1 

for ^ + Q + Sr — 7r/4 + 7tv, v — 0, 1, ... and its minimum 

^ J = 2ka 6 I (18) 

* n / min 

for | + Vt + 5 R = 3vr/4 + vrw, v = 0, 1, ... . 

Let us now turn to the opposite limit of strong absorption, characterized by Si ^> 1. In 
this case we get 

a = a (l — i), (19) 

r e = 0, (20) 

(Tel = 8jca 2 , (21) 

(Tin = ^j- (22) 

The above cross-sections (eqs. I2TH221) are determined by the van der Waals tail of the 
adiabatic potential only. Remarkably, they include no information about the short- and 
middle-range parts of the atom-antiatom interaction. All particles that penetrate to these 
short distances are lost, so the only information available to the observer is due to the 



so called quantum reflection [18|, ll9j from the asymptotic tail of the potential, which is 
determined by C§ alone. A remarkable feature of the strong absorption limit is that the 
effective range ([9]) is exactly zero. 

In Fig. [1] we plot the elastic cross-section as a function of real part of short-range phase- 
shift 5r for two values of the "inelasticity parameter" 5j = 0.3 and 6j = 0.7. Pronounced 
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- 5=0.7(modcl III) 




FIG. 1: Elastic cross-section for H—H collisions as a function of phase Sr, calculated at E — > 
according to a e i = 4-/r|a| 2 with a given by eq. ([7]). 




FIG. 2: Annihilation cross-section for H — H collisions as a function of phase 5r, calculated at the 
energy E=10 -6 a.u. according to <7j n = A-Klm{a)/k with a given by eq. ([7]). 

oscillations can be seen in the "weak absorption" case, which become much more smooth 
with increasing inelasticity parameter (in the limit 5j ^> 1 the cross-sections are constant). 
The same tendency is apparent in the inelastic cross-section (Fig. [2]) and in the ratio of the 
elastic to the annihilation cross-section (Fig. [3]). 




S=0.3(modell) 
5=0.7(modcl III) 



FIG. 3: Ratio of the elastic to the annihilation cross-section as a function of phase 5r, calculated 
at the energy E=10 -6 a.u. with a given by eq. (J7J). 

III. PHYSICAL IMPLICATIONS 

In this section we discuss the physical implications of the sensitive dependence of the 
scattering length on the short-range phase shift. We study the evolution of the scattering 
length as function of the phase-shift 5r. This can be considered as an universal tool for 
studying the contribution of various physical effects to scattering observables. 



1. Limits on the accuracy of the model His — His calculations 



Since the complete ab initio treatment of the hydrogen- ant ihydrogen scattering is a diffi- 
cult problem it is important to be able to estimate the accuracy of the approximate calcula- 
tions. Several model calculations have been performed to get elastic and annihilation cross- 
sections in the case of H lS - H 13 scattering [3, \i, U U, Id, bo|, M, U2, bd, hi ha, ha, IzZI, bsj- 
In these models mostly the same adiabatic potential V a d of H is — His was used, applicable 
in the distance range r > R s ~ \tb [lJ]- At the same time very different approximations 
were made regarding the rearrangement and strong forces. They include calculations of the 
elastic cross-section neglecting both rearrangement and strong force 7( , the account of rear- 



6, 
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29] , calculation of the annihilation 



rangement in the " distorted- wave" approximation 

cross-section due to strong force only [2jJ, [23(, non-perturbative calculations based on the 

approximative optical potential p, IS] and other. 



In terms of the formalism developed in the previous section the above mentioned models 
differ by the short-range phase S = Sr + Si, i.e. the phase accumulated up to R ~ R s . This 
phase incorporates the effects of inelasticity due to strong forces and rearrangement. In the 
optical model calculations by Voronin and Carbonel [5j (further referred to as model I) it 
was found that a = 5.2 — il.8 a.u. (without account of strong force). From this value of 
the scattering length the corresponding value of the short-range phase-shift can be uniquely 
deduced (in the sense that the solution of the Schrodinger equation with given adiabatic 
potential V a d and with boundary condition flTJ results in unique correspondence between S 
and a, as can be seen from eq. (Cj)). The value of the short range phase-shift for model 



I turns out to be S = 0.696 + i0.305. The value presented in [22j by Armour, Liu and 
Vigier (further referred to as model II) is a = 8.2 — z2.8 a.u. and the short-range phase-shift 
deduced from this value of the scattering length is S = 0.136 + iOAl. The optical model 
calculation |8j by Zygelman, Saenz, Froelich and Jonsell (further referred to as model III) 
gives a = 5.6 — i3.7 a.u. (without strong force). The corresponding short-range phase-shift 
turns out to be S = 0.638 + z0.715. Finally, the calculation of elastic scattering neglecting 



both rearrangement and strong force effects (further referred to as model IV, [29]) results 
in the scattering length a = 7.69 a.u.. The corresponding real phase-shift, deduced from 
this value of the scattering length is S = 0.287. All these values are tabulated in Table [B 
Inclusion of the inelasticity is thus equivalent to introduction of the short-phase corrections 
to model IV. Whereas it was possible to separately calculate the complex phase due to direct 



annihilation 
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30J], calculation of such phase correction due to the rearrangement is still 
incomplete since it is based only on the imaginary component of the optical potential 8(. 
Hence the inaccuracy of the real phase Sr remains an important source of the inaccuracy of 
the scattering length in H — H scattering. 

From the eqs. (I15II16P it is clear that the smaller is the inelasticity parameter Si, the 
more sensitive are the model results to the uncertainty in the real part of the H — H 
interaction at distances r < R s . We notice that for a fixed value of inelasticity parameter Si 
the relation between imaginary and real part of the scattering length can not be arbitrary 
but is determined by eq. ([7j). To illustrate this statement we plot the possible values of the 
scattering length a in the complex plane of a as a function of Sr (which changes from —tc/2 
to 7r/2) for fixed value of Si. The possible values of a for each fixed value of Si form closed 
curves. The three curves of the scattering length correspond to three different inelasticity 

10 




5=0.3(model I) 

— S=0.4(raodel II) 
5=0.7(model III) 
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FIG. 4: Complex scattering length a as function of the short-range phase-shift 5r, for the fixed 
value of 5j. The square indicates the complex value of the scattering length obtained in model 
I, circle - in model II, triangle - in model III. The value of a ~ 5 — z5 a.u. marked by a star 
corresponds to the limiting case of strong absorption 8i 3> 1. 



parameters 5j deduced from models I-III. One can see that the smaller is 5j the larger is 
the possible variation of a with Sr. In the limit of zero absorption (5j — > 0) the variation 
becomes infinite and the closed curve degenerates into the whole real axis of the complex 
plane (model IV), while in the opposite limit of strong absorbtion the curve shrinks into a 
point a = ao(l — i) — 5 — z'5 a.u.. Interestingly, the mean value of the scattering length 
a (averaged over Sr) is the same for all curves, i.e. it turns out to be independent on the 
inelasticity parameter 5j and is equal to a — a (l — i). 

The corresponding variations of the cross-sections are shown on Fig. [1] and Fig. [2] for 
models I and III (the corresponding inelasticity parameters are Si = 0.3 and Si = 0.7, 
respectively). One can see that the ratio of maximum to minimum possible values of elastic 
cross-section turns out to be cr^ a;E /(T™ n ~ 46 for Si = 0.3. The same value for Si = 0.7 
is much less and turns out to be cr™ ax /a™ m ~ 4. In view of this result it is clear that the 
knowledge of the absorptive phase-shift Si is very essential and it is sufficient for establishing 
the limits of accuracy in calculations of both inelastic and elastic cross-sections. 

We can perform further analysis and estimate how the uncertainty in the H — H interac- 
tion at distances r < R s would be reflected in the uncertainty of the scattering observables. 
Such an estimation can not be model independent. We will use the following simple assump- 
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tions. One can expect 2l|, |22|, |28| that the leading contribution to the mentioned interaction 



at distances R < R s is given by the local interaction of the form (expressed in a.u.): 

V loc (r) = -1/r + 0.75. (23) 

This potential continuously matches with V a d{r) at r m R s a.u. Apart from such local 
potential there is an additional interaction in the H — H channel, which results from coupling 
to the decay channels. We will model the effects of coupling to the Pn and Ps channels, 
responsible for rearrangement, by the local complex potential of the form: 

V ovt = ^^lexp(-2r). (24) 

The above simple form of "optical" potential is only used to mimic the effect of coupling 
to the decay channels and can not be treated as true form of effective short-range H — H 
interaction - but it captures its localization. The ab initio calculations of complex nonlocal 
optica! potential (so far restricted by app—ons) can b e round in fl fa Tirese 
optical potentials are constructed by projecting the Hamiltonian on the subspace of open 
Pn and Ps channels and therefore their characteristic range is equal to the "size" of the Pn 
states with energetically highest possible principal quantum numbers (N < 24). Because of 
that such ab initio optical potentials are localized below R s and this property is reflected in 
the r dependence of the simple model potential V opt {r). 

To calculate the short-range phase-shift we will use the semiclassical approximation for 
the phase: 

5^ f S ^/-2M[V loc (r) + V opt (r)}dr - vr/4. (25) 

The term 7r/4 is the quantum correction to account for the correct behavior of the coulomb 
wave-function at small distances r < 1/M, where semiclassical approximation is no longer 
valid. Expecting |K pi (r)| <^ |VJ oc (r)| we get from (j25|) for 5: 

5 = 5° + 5 opt , (26) 

rRs 

5° = / y/-2MV loc (r)dr, (27) 

Jo 

S opt = _^ f Rs Vopt(r) dr (2g) 

Jo ^-2V loc {r) 

In the above expression 5° is a phase-shift produced by the local potential VJ oc alone, while 
5 opt is the variation of that phase-shift due to the presence of the optical potential. Using 
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the adopted expressions for VJ oc and V opt we come to the following numbers: 5° = 72.45, 
6 o P t = _( v _i w )6.38. 

The corresponding variation of the scattering length can be obtained from eq. (1141) : 

5 opt 

6a = ~ a ° sin 2 (7r/8 + n + (yo + ^)- (29) 

From the above numbers we can conclude, that the value of inelasticity parameter 5i consis- 
tent with models I-III corresponds to the value of the imaginary part of the optical potential 
w which lies in-between 0.05 and 0.1 a.u. 

On Fig. [5] and Fig. [6] we plot the elastic and inelastic cross- sect ions as function of the 
real part of the model optical potential V opt . Even a small and confined real part produces 
large variations of the cross sections. For the optical potential given by eq. (124")) . the 
mean variation of the ratio v opt within the interval [0, R s ] is on the order of 1%, yet it 
produces more than 100% change of elastic and inelastic scattering-cross sections. This 
result makes clear the relative importance of the real part of the effective H — H short- 
range interaction and puts obvious restrictions on the accuracy of calculations. It means 
in particular that the non-perturbative account of additional decay channels could result in 
significant changes of the elastic and inelastic cross-sections. We will show in the following 
that such an "instability" of model calculations (for relatively small values of inelasticity 
parameter Si) is due to the existence of the spectrum of narrow near-threshold states in the 
H — H channel. 

2. Sensitivity to the nuclear interaction 

An important consequence of the weak absorption in His ~ His interaction is the possi- 
bility of observing the nuclear effects on the molecular scale. The effect of the nuclear forces 
was reported in ref. [5j, |2l|, 
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271 ] for the case of hydrogen-antihydrogen scattering and 
in ref. |30| for antihydrogen-helium scattering. In calculations 2jJ, |22J no rearrangement 
effects were taken into account. 

The strong force effects can be studied by inclusion of the complex short-range nuclear 
potential in the nonrelativistic Hamiltonian describing the H H system. This potential ac- 
counts for the nuclear interaction of proton and antiproton on the scale of few fm (1.88 • 10 -5 
a.u.) and is chosen to reproduce the nuclear pp (spin dependent) scattering length. Incorpo- 
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FIG. 5: Elastic cross-section as a function of the real part v of the optical potential (I24|) . calculated 
at the energy E=10 -6 a.u. 
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FIG. 6: Inelastic cross-section as function of the real part v of the optical potential (I24j) . calculated 
at the energy E=10~ 6 a.u. 

ration of the short-range complex potential in the model calculations results in modification 
of 5r and 5j. Strictly speaking, the effect of strong forces is not restricted to "direct" an- 
nihilation and scattering on the nuclear potential in the H — H channel. It appears also 



through modification of the final states and energies of Protonium. Neglecting 



ment these "second order" effects, one can use the well-known expression 31 



or the mo- 



32| for the 



nuclear phase-shifts in the presence of Coulomb interaction to extract directly the phase-shift 
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of interest: 

k C(rj) cot 5 sc + 2r]h(i]) = -l/a sc 

where a sc is the Coulomb corrected nuclear scattering length, rj = —M/k (k is the cm. 
proton-antiproton momentum) and functions C(rj), h(rj) are given by 

C{V) = 77; — v T, 

exp{2irri) — 1 

h{rj) = l/2[W(-iri) + V(irj) - hi(r/ 2 )]. 

In the above, \I/(x) = T'(x)/T(x), where F(x) is a standard gamma-function. Neglecting 
the leptonic energy in comparison with the value of Coulomb pp potential at the distance 
around 1 fm, we can put k — > and get: 

(L w -2vrMa., r . 



The values of the strong-force scatering lenghth a sc in the Kohno-Weise model 33|, |34J are: 



d sc (S = 0) = (1.07-*1.45)-10~ 5 a.u., (30) 

a sc (S = 1) = (1.68 - zl.06) ■ 10" 5 a.u. (31) 

while the corresponding nuclear phase-shifts are 

8 SC (S = 0) = -0.06 + i0.08, (32) 

5 SC (S = 1) = -0.1 + i0.06 (33) 

where S = 0, 1 are the values of total nuclear spin. 

The models I- IV include no information about strong forces. However we will show 
that, once the value of scattering length without any account of strong forces is known, one 
can also obtain the corrected value of this scattering length that includes the presence of 
strong forces. In fact, adding the nuclear phase-shifts (1321) and (1331) to the already discussed 
phase-shifts deduced from the rearrangement calculation of models I-IV 

5 -> 5 + S sc , 

one gets the values of the scattering lengths and cross-sections, corresponding to simulta- 
neous account of rearrangement and strong force effects. The above mentioned addition of 
phase-shifts can be done due to the fact that contributions to the phase from each of the 
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Model I 


Model II 


Model III 


Model IV 


a, a.u. 


5.2-11.8 


8.242.8 


5.643.7 


7.740 


a sc (S = 0),a.u. 


5.642.2 


8.343.4 


5.844.0 


8.440.4 


a sc {S = 1), a.u. 


5.842.1 


8.643.4 


6.043.9 


8.140.6 


6 


0.696+10.305 


0.136+i0.41 


0.638+i0.715 


0.287+iO 


a e i, a.u. 2 


380 


943 


566 


745 


a s J(S = 0), a.u. 2 


448 


1018 


623 


889 


a s J(S = l), a.u. 2 


472 


1077 


640 


829 


VEa in , a.u. 


0.53 


0.8 


1.08 





^£< c (S = 0),a.u. 


0.65 


1.0 


1.17 


0.13 


^£< c (S = l),a.u. 


0.62 


1.0 


1.15 


0.18 



TABLE I: Values of the short-range scattering lengths, phase-shifts and cross-sections (in a.u. 2 ) 
calculated at the energy E = 10~ 6 a.u. with the parameters deduced from different models I-IV 
for H — H interaction. The superscript sc corresponds to the account of the strong force in the 
presence of Coulomb field. The results were obtained by numerical solution of the one channel 
Schrodinger equation subject to the boundary condition given by eq. Q]with the phase 5 from row 
5 incremented by the strong force contribution given by eqs. 



effects (strong force and rearrangement) come from very different distances. This is true 
until the effect of strong force on the final Pn states is not taken into account. 

The results of our calculations for the elastic and inelastic scattering cross-sections with 
and without strong force, based on the parameters deduced from models I-IV are collected in 
Table [H We observe that the simultaneous account of rearrangement and strong-force effects 
significantly differs from the results where only one of these effects is taken into account. 

The modification of cross-sections due to the strong force is significant when compared to 
the ratio of nuclear and atomic scales, which is on the order of 10~ 5 . One can also see that 
this modification is model dependent. This fact is clear from our previous considerations of 
the cross-sections as function of the the short-range phase-shift. In particular, the change 
of the elastic cross-section due to the strong force is around 24% in model I, while it is 
13% in model III. These models differ by the inelasticity parameter 5j which is 2.3 times 
larger in model III as compared to model I. Consequently model III is less sensitive to any 



1(3 



variation of phase, in particular the phase induced by strong interaction. The physical reason 
for the lack of sensitivity to the details of the short-range interaction with the increase of 
inelasticity parameter is clear. Indeed, the amplitude of the wave reflected back into the 
H — H channel that is generated at such small distances and "carries back" information 
about the nuclear forces is, in the case of strong absorption, exponentially small. Hence in 
the case of strong absorption the scattering length only weakly depends on the details of the 
short range interaction. 

Let us mention here that the effect of nuclear forces on H — H interaction would be 
strongly enhanced in the vicinity of certain value of the real part of the short-range phase. 
According to (1141) the largest variation of the scattering length with small variation of the 
phase 5r is achieved when the following condition is fulfilled 

tt/8 + ft + 5 R = 7cv (34) 

where v is an integer number. As shown later in sec. IHI 41 this condition corresponds to the 
appearance of the loosely bound H H state. In particular, for the inelasticity parameter of 
model I (5j = 0.3) the modification of the elastic cross-section due to the strong force in the 
near-resonance situation (eq. [M]) would be from a e i = 3800 a.u. 2 without strong force to 
a s e f(S = 1) = 2015 a.u. 2 and a s J(S = 0) = 2249 with strong force, i.e. around 50%. 

3. Isotope effect 

The oscillatory behavior of the cross sections gives rise to the strong isotope effect. Simple 
changing of the reduced mass alters the phases 5 and Q in eq. ffl4l . The mass dependence of 
the phases can be easily established in the WKB approximation (l25l) . from which it follows 
that the semiclassical phase is scaled as the square root of the reduced mass M. The change 
of the reduced mass from M p /2 to 2/3M p (where M p is the proton mass) corresponds to 
the replacement of hydrogen by deuterium and results in multiplication of the semiclassical 



phase by the factor a/4/3 



n d = a/473^, 

5 d = y^/3(6 - vr/4) + tt/4. 
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Here the subscript d refers to deuterium. In the above formula we took into account the 
Coulomb quantum correction — 7r/4. 

The numerical calculation of the scattering length while neglecting rearrangement and 
strong force (i.e. using potential VJ oc (r) as the D — H interaction for r < R s ) gives 

a DR = -41.07 a.u. 

to be compared with the corresponding scattering length (without rearrangement and strong 
force) in the H — H case 

a HR = 7.69 a.u. 

The account of rearrangement effects dramatically changes the situation. 

The model calculations by Voronin and Carbonell [5j show that the scattering length 
changes from a = 5.2.1 — zl.8 a.u. for H — H to 15.0 — ill. 6 a.u. for the H — D case. In 
view of the demonstrated high sensitivity of the results to the short range interaction, very 
accurate calculations are required to obtain the reliable values of the cross-sections. However 
the prediction of a strong isotope effect has a general character and is a consequence of "weak 
absorption" in the H — H interaction. Note that in the opposite case of "strong absorption" 
the expression (fT9l predicts only weak monotonous dependence of the scattering length on 



the reduced mass via a ~ v / 2MC 6 . 

4- The near-threshold quasimolecular states of H — H 

The oscillating behavior of the scattering cross-section as a function of the phase 8 has a 
clear physical meaning. The long range character of the attractive atom-antiatom interaction 
combined with the weak absorption supports the existence of long-lived highly excited H — H 



states 



2o[ 



26] . They manifest themselves as near-threshold singularities of the S'-matrix 
and thus strongly affect the low energy scattering. By changing the phase through varying 
the strength of interaction, changing the reduced mass of the system, including additional 
decay channels, etc., the position of such states can be tuned to be at the threshold, inducing 
a rapid increase of the cross-sections as a function of the tuning parameter. The S'-matrix 
of the finite range potential can be expressed through the position of all its poles via the 
following expression 35J: 

S = ] T -exp(—2ika) (35) 

k 
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where Zk is the pole of S-matrix and a is independent on k. 

If there is a S-matrix pole z, situated much closer to the threshold than other poles, we 
can approximate the infinite product by one term only: 

S(k) « ^—exp(-2ika). (36) 

z — k 

In the limit of small k the expression (j36p becomes S(k) = 1 + 2k/ z — 2ika, from where one 
gets the connection between the scattering length and z 

a = i/z + a. (37) 

When there are two (or more) poles at similar distance from the threshold (bound and virtual 
states) one should take into account the contribution of the number of equidistant poles of 



the S'-matrix so that t 
It was shown in 12 



re expression for the scattering lengths becomes a = i J2k l/ z k + a - 



36J that the above expression is applicable in case of scattering on 
potentials with 1/r 6 tail. In such a case the constant a in the one-pole expansion is a = ao, 
with a given by eq. (|Sj). Therefore the position of the pole in the one-pole expansion of the 
S'-matrix is given by 

z = i/(a — a ). 

This relation is valid in the case when the distance between poles is much bigger then z. 
Substituting the expression for the scattering length given by eq. (JTj) we get 

z= — tan(<5 + fi + 7r/8). (38) 

ao 

Taking into account that the near-threshold states are possible only in case of weak absorp- 
tion ((5/ < 1) we obtain 

z = —tan{5 B + Q + 7i/8) ,._ l — . (39) 

a a cos 2 {5 R + Q + 7f/8) v ' 

The above equation is valid only for the expansion with respect to the single nearest- 
threshold pole, which gives the dominant contribution to the scattering length. 

From the above equation it is clear that due to the inelasticity of scattering characterized 
by 5j, the near-threshold S-matrix poles are shifted to the left from the imaginary axis of 
the complex /c-plane. These shifted poles in the 2-nd quadrant of the complex fc-plane can 
be interpreted as localized (bound) states with an inelastic width, while the poles in the 3-rd 
quadrant are the virtual states or resonances with inelastic width. (The distinction between 
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resonances and virtual states can be obtained by a continuous decrease of inelasticity. In the 
limit of zero inelasticity the virtual states come to the negative imaginary axis of the /c-plane, 
while the resonances are manifested as pairs of poles in the 3-rd and 4-th quadrants.) 

The condition for the appearance of the new state now means that the pole comes from 
the third to the second quadrant of the complex fc-plane. At the moment when the pole 
of the S-matrix crosses the real (negative) axis of the complex k-plane one has Imz = 0. 
According to eq. (139]) this is equivalent to condition (|34j) : 

5 R + VI + tt/8 = nv 

where v is an integer number. Numerical integration of the adiabatic potential V a d gives the 

WKB phase Q = 19.383 so that the numerical value of 5r (by modulus it) corresponding to 

the resonance is 5 R = —0.926. 

We have shown that the previously introduced "resonance" condition ( 1341) is indeed the 

condition for the appearance of a new state at the threshold. At the moment of appearance 

of the new state the scattering length is: 

/ cosh(5/)\ , A . 

a res = o 1 - i . r r? x 40 

V sinh(dj) / 

and its imaginary part attains its maximum value. In case of weak absorption (5j <C 1) it 

can be written as 

Ima res = --. (41) 

At the same time, the variation Aa due to a small change of the real part of the phase-shift 

A5r attains its maximum value: 

__ 2a A5 R 

cosh(5/) 
Hence the effect of additional interactions (e.g. the account of additional decay channels, 

strong force, etc.) in the near-resonant situation will result in a fast change of the real part 

of the scattering length. The point at which the S'-matrix pole crosses the real axis of the 

complex fc-plane can be obtained by comparing formulas (|39|) and (|4T|) : 

%res 

a 
It is determined entirely by the inelasticity parameter. The larger is the inelasticity param- 
eter, the larger is the shift of the S'-matrix pole to the left from the origin of the complex k 
-plane, and the smaller is its influence on the scattering cross-sections. 
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The energy of the near-threshold state is: 

z 2 tan 2 (<J fl + VL + vr/8) i<5 7 tan(5 R + Q + ir/8) 

b ~2M~ 2Ma\ Ma 2 cos 2 {5 R + Q + 7r/8)' ^ ' 

If the pole z is very close to the real axis of the complex /c— plane, so that Sr+Q+tt/8 = nv+tp 
and ip <C 61, the expression (H2l) should be modified: 

Eo " 2Ma 2 ■ (43) 

Here ip measures the de-tuning of the phase from the exact resonance. The state at the 
threshold is extremely extended spatially and has a very small width: 

r 2yfr 

2Ma 2 ' 

The detailed study of the near-threshold states in 1/R e potential can be found in [l7|, |36j. 
In the above cited papers the higher order expansion in ao/(a — ao) for the energy of near- 
threshold state is given: 

1 / 0.918a 0.947a 2 , 

-^0 = -77TT7 ^7 1 + 



2M(a — ao) 2 \ a — ao (a — ao) 
Taking into account (j7j) we get: 

E = ~ tan ^77 2 + ( X + 0.918 tan(5 + Q + vr/8) - 0.947 tan 2 (5 + Q + vr/8)) (44) 

We observe that once the scattering length a is fixed by a certain model, the spectrum 
of near-threshold quasi-molecular states is also fixed (up to the total number of H — H 
states with given hadron angular momentum). Indeed, the value of the scattering length 
uniquely determines the short-range complex phase-shift 5 and thus the boundary condition 
(Crj, while the solution of the Schrodinger equation with given adiabatic potential and fixed 
boundary condition ([1]) gives the spectrum (up to the number of states). The validity of this 
statement is restricted only by the approximation that our boundary condition ([1]) is energy 
independent, i.e. it is valid for energies E <^ e m ~ 0.05 a.u.. We calculate the positions 
of the bound states nearest to the threshold, as they are fixed by the scattering lengths of 
various models. These values are collected in Table [TT1 We also present the position of the 
same states when the spin dependent strong force is taken into account. One can see that 
the shift of the states due to the account of nuclear forces is model dependent and is of order 
10~ 5 a.u. for the first loosely bound state. The splitting AE SC of the triplet and singlet spin 
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Model I 


Model II 


Model III 


Model IV 


E, 10" 5 a.u. 


— 10.1 — z5.6 


-2.2 - z3.2 


-3.8 - il3.1 


-4.9 - i0 


E SC (S = 0)10~ 5 a.u. 


-8.4-J6.7 


-1.4 - z3.3 


-0.7-il2.8 


-4.3-t0.8 


E SC (S = 1)10~ 5 a.u. 


-7.8 — J6.1 


-1.3-i2.9 


-0.8 -ill. 5 


-3.9-i0.6 



TABLE II: Energy of the state nearest to the threshold. The first row collects the values without 
account of the strong force, the second and the third row collects the values with account of the 
strong force in spin states 5 = and S = 1 correspondingly. The results were obtained by numerical 
solution of the one channel Schrodinger equation subject to the boundary condition given by eq. 
[1] with the phase including the strong force contribution, together with the square-integrability 
requirement. 

states is on the order of 10~ 6 a.u., which is one order less than the inelastic width due to 
rearrangement. 

It is useful to see, how the position of the pole that is nearest to the threshold changes 
with the short-range phase-shift S. Such an evolution of the poles in models I-III is presented 
on Figd Each of the shown trajectories corresponds to a certain fixed value of inelasticity 
parameter Si. The evolution of the S- matrix pole along the trajectory can be understood as 
a gradual " switching on" the real part of effective interaction, responsible for rearrangement 
(while its imaginary part is kept fixed). At certain value of the phase Sr the pole crosses 
the real axis of the k-plane, and the new vibrational state (quasi-bound) appears in the 
spectrum of H — H system. With increasing the real part of effective interaction (and 
correspondingly the phase Sr) this pole goes far away from the threshold, becoming deeper 
bound state. At the same time a new pole with the negative imaginary momentum (virtual 
state) is approaching the real axis. This is reflected in oscillatory behavior of the elastic and 
inelastic cross-sections. 

We conclude by emphasizing that the resonant dependence of the cross-sections on the 
interaction parameters (such as real and imaginary short-range phase-shift) is of purely 
quantum mechanical origin and occurs only in the case of weak absorption. This resonance 
behavior reflects the existence of the spectrum of narrow, near-threshold, quasi-molecular 
states of H — H. 
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Re k, (a.u 



FIG. 7: Evolution of the poles of the S'-matrix in the complex £;-plane as function of the real part 
of the short-range phase-shift, —tt/2 < 5r < n/2. Solid line corresponds to the fixed inelasticity 
parameter 5i = 0.3 (model I), dash-dotted line - inelasticity parameter 5i = 0.41 (model II) and 
dashd line - inelasticity parameter 5j = 0.7 (model III). An empty square marks the virtual state, a 
star marks the point where the S-matrix pole crosses the real axis of the complex k-plane, a circle 
- the position of the bound state corresponding to the "true" value of the phase-shift 6r, deduced 
for each model from the corresponding scattering length. In calculations, the diverging character 
of the virtual states was tempered using the technique of exterior complex scaling [371 ] . The arrow 
shows the direction of the pole evolution with increasing 5r. 



IV. CONCLUSIONS 

We have presented an overview of the H — H collisions based on the analytical treatment 
of the asymptotic potential tail in conjunction with the parametrization of the contribution 
from the inner part of the interaction. We show analytically that for a realistic value of the 
inelasticity parameter 5j (the imaginary part of the short-range scattering phase induced by 
the rearrangement and the strong-force interaction) a rich spectrum of narrow near-threshold 
HH states exists. These states, corresponding to the near-threshold singularities of the S'- 
matrix, determine the behavior of the elastic and inelastic cross-sections. The latter exhibit 
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high sensitivity to even small perturbation of the H — H interaction or its approximate 
treatment. We predict important physical effects such as a strong isotope effect (i.e. the 
resonance-like dependence of elastic and inelastic cross-sections on the reduced mass of the 
system) and demonstrate a significant dependence of both elastic and inelastic H — H cross- 
sections on the nuclear potential. 
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